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ABSTRACT 


We used a novel, overdamped oscillator aboard the Space Shuttle to measure the viscosity 
77 of xenon near its critical density p c and temperature T c . In microgravity, useful data were 
obtained within 0.1 mK of T c , corresponding to a reduced temperature t = (T - T c ) /T c = 
3 x 10 . The data extend two decades closer to T c than the best ground measurements, and 

they directly reveal the expected power-law behavior 7 ] oc t~ wz *. Here v is the correlation 
length exponent, and our result for the small viscosity exponent is z, = 0.0690 ±0.0006. (All 
uncertainties are one standard uncertainty.) Our value for z ^ depends only weakly on the 
form of the viscosity crossover function, and it agrees with the value 0.067 ± 0.002 obtained 
from a recent two-loop perturbation expansion [H. Hao, R.A. Ferrell, and J.K. Bhattacharjee, 
Preprint (1997)]. The measurements spanned the frequency range 2 Hz < / < 12 Hz and 
revealed viscoelasticity when t < 10~ 5 , further from T c than predicted. The viscoelasticity 
scales as Afr , where r is the fluctuation-decay time. The fitted value of the viscoelastic time- 
scale parameter A is 2.0 ± 0.3 times the result of a one-loop perturbation calculation. Near 
T c , the xenon s calculated time constant for thermal diffusion exceeded days. Nevertheless, 
the viscosity results were independent of the xenon’s temperature history, indicating that 
the density was kept near p c by judicious choices of the temperature vs. time program. 
Deliberately bad choices led to large density inhomogeneities. At t > 10 -5 , the xenon 
approached equilibrium much faster than expected, suggesting that convection driven by 
microgravity and by electric fields slowly stirred the sample. 

PACS index categories: 51.20. +d; 05.70Jk; 64.60.Ht; 83.50.Fc; 83.85.Jn 
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1. Introduction 


As the liquid- vapor critical point is approached, the shear viscosity rj (£) measured in the limit 
of zero frequency diverges as £**>, where £ is the correlation length, which itself diverges on the 
critical isochore as £ - “. (Here, t = (T — T c ), T c is the critical temperature, and u = 0.630.) 
The viscosity exponent z n is central to the theory of dynamic critical phenomena [1]; the 
quantity 3 + z v is called the ” dynamic critical exponent” and it characterizes transport 
phenomena in all near-critical fluids with a scalar order parameter. Because the exponent 
z,, = 0.069 is so small, it is very difficult to measure accurately on Earth in a pure fluid 
such as xenon. Far from the critical point, the separation of the viscosity’s divergent critical 
contribution from the analytic background contribution is uncertain because the critical 
contribution is a small fraction of the total and because the separation depends sensitively 
upon the “crossover function” H (£), which is known only approximately. Close to the critical 
point, the very compressible xenon stratifies in the Earth’s gravity and the divergence of 
the viscosity is blunted in a manner that depends upon the height of the viscometer. As 
illustrated by Figure 1, stratification visibly blunted the divergence of the viscosity of xenon 
near £ min = 3 x 10 -5 , even when the measurements were made using a high-Q oscillator only 
0.7 mm high. The arguments of Ref. [2] show that £ min scales with the height h and the 
acceleration of gravity g as (g/i) 0 64 . Thus, a significant reduction of £ min requires either a 
viscometer with h << 0.7 mm on Earth or the use of microgravity. 

Here we report results obtained with a novel viscometer that we integrated into the 
“Critical Viscosity of Xenon” (CVX) experiment package and operated aboard the Space 
Shuttle Mission STS-85. In the Shuttle’s microgravity environment (approximately 1 x 10 -5 
ms -2 ), CVX obtained useful viscosity data at £ inin = 3x 10 -7 , two decades in t closer to 
T c than the best ground-based measurements. Indeed, the present data are useful closer 
to T c than the data from any prior studies of liquid-vapor critical points, including those 
conducted in microgravity ( e.g . [3, 4, 5) [6, 7]). The CVX data yielded a more accurate 
value of z n than previous work, and it also produced the first accurate measurements of the 
frequency dependence of viscosity near any critical point. 
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The CVX viscometer was a torsion oscillator that was driven at frequencies in the r an ge 
2 Hz < / < 12 Hz. Both the in-phase and the quadrature components of the oscillator’s 
response to the driving torque were measured. This permitted us to separately deduce the 
real and the imaginary components of the viscosity, Re (77) and Im (77). Both components are 
required to fully test theory. The imaginary component is the elastic part of the response 
to shear stress that is ordinarily studied either at much higher frequencies or in complicated 
fluids such as polymer melts. In near-critical fluids, Im (77) ^ 0 when fr> 1, where / is the 
frequency of the measurement and r is the fluctuation-decay time 


= P$- = t 0 r^>. 


^sT c 


( 1 . 1 ) 


where kg is Boltzmann’s constant. At the very low frequencies used by CVX, the viscoelastic 
behavior of xenon was evident when t = 10 -5 and it was dominant when t — 10 -6 . (See 
Figure 2.) 

The theoretical model used to analyze the CVX data combined the scaling function S (z) 
for near-critical viscoelasticity from Ref. [8], the crossover function H (£) from Ref. [9], and 
the background viscosity 770 from Appendix C of Ref. [10] to obtain the prediction 


V ((,/) = Vo exp [z,//(0][S(Az)]"^ 3+ ^. (1.2) 

In Eq.(1.2), the scaled frequency is defined by z = -nr/r, and we introduced the parameter 
A into the argument of the scaling function S (z) to obtain agreement with our data. Accu- 
rate measurements of the scaling function require accurate measurements of 77 (£, /) in the 
region where is large and where z is varied through a large range from z « 1 to 2 » 1. 
Measurements prior to CVX were unable to achieve these conditions [11, 12, 13, 14]. 

The CVX viscosity data at 2 Hz, 3 Hz, 5 Hz, 8 Hz, and 12 Hz were used to determine 
five parameters. Two are the “universal” parameters z tf and A\ two are the wave vectors q c 
and q D that occur in the crossover function for xenon; and one is the value of T c on CVX’s 
temperature scale. When all of the data within the range 10~ G < t < 10“ 4 were fitted for 
these parameters, we obtained z v = 0.0690±0.0006, q c = 0.051±0.007, and q D = 0.16±0.05. 
(The uncertainties indicated throughout this report, are one standard uncertainty and they 
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As discussed in Section 5.3, the CVX data determined the product At q = 2.31 ± 0.06 ps 
with a relative uncertainty of 0.03. The relative uncertainty of A is larger (0.15) because of 
the uncertainty of the correlation length amplitude £ 0 that propagated into r 0 and then into 
Atq. The uncertainty A could be reduced by a factor of five if the uncertainty of £o were 
reduced, perhaps by additional measurements on Earth. 

The present research was preceded by viscosity studies on Earth of xenon and carbon 
dioxide near their liquid-vapor critical points [10] and binary liquid mixtures near their 
consolute points [14]. The latter belong to the same dynamic universality class as pure fluids 
but are less influenced by gravity. We now identify the ways that CVX complements the 
previous research. 

For tests of theory, i.e. Eq.(1.2), a simple pure fluid such as xenon has three advantages 
in comparison with binary liquid mixtures near consolute points. First, the dependences 
of the noncritical viscosity tj 0 (p, T) on density p and temperature T are weak compared 
with the dependencies of r) 0 ( x , T ) on mole fraction x and temperature for binary liquids. 
For pure fluids, the temperature dependence is indistinguishable from the well-understood 
dilute gas behavior 770 (0, T) [19]. This means that r) 0 (p,T) has a temperature dependence 
that is weaker than for binary liquids, and its use in fits to Eq.(1.2) adds no free parameters. 
Second, the fluid-dependent parameter qc that appears in the crossover function H (£) [9] is 
known for xenon; it has not been determined for any binary liquid. (The other parameter 
in H (£) is the effective cutoff wave vector qo that is used in a mode coupling integral over 
momentum space. It is a free parameter related to the amplitude of the viscosity divergence.) 
Third, current technology allows the conditions of low frequency and shear rate to be met 
more easily with pure fluids than with liquid mixtures because the decay time r is at least 
100 times faster in pure fluids than in mixtures at the same reduced temperature. 

Major divisions of the remainder of this manuscript are titled: Apparatus, Temperature 
timeline, Data reduction, and Results. Appendices deal with: Tabulated viscosity data, 
Electrostriction, Frequency-dependent scaling function S ( 2 ), and Estimation of qc- 
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2. Apparatus 

The flight apparatus consisted of the cell holding the xenon sample and the oscillator, the 
surrounding thermostat, and the electronics, all contained in two flight cannisters. Here we 
provide an overview. Additional details can be found in References [17] and [18]. 

2.1. Oscillator and sample cell 

The heart of the CVX viscometer, shown in Figure 5, was an oscillator contained in a thick- 
walled copper sample cell. The oscillator was an 8 x 19 mm rectangle of screen that was cut 
out of a larger piece of nickel screen. The screen consisted of 0.03-mm wide wires formed 
by electrodeposition (Buckbee-Mears [20]) and spaced 0.85 mm apart. When the screen was 
cut, two wires were left extending from the edges of the rectangle. They served as torsion 
wires and were soldered to a stiff yoke that was centered between four electrodes parallel to 
the screen. We chose the oscillator’s dimensions so that it could be assembled by hand and 
so that it would be sensitive to changes of viscosity at low frequency. The oscillator’s mass 
was approximately 1 mg. 

Four electrodes were used to apply torque to the oscillator. To do so, diagonally opposite 
pairs of electrodes were charged to different voltages while the screen remained grounded. 
The electrodes were 10 x 12 mm rectangles of 0.13 mm thick brass sheet soldered to 1 mm 
diameter wires that acted as both mechanical supports and as electrical leads. The electrodes 
lay in two planes located approximately 4 mm above and below the plane of the oscillator. 
The screen s yoke and the four electrodes were supported by feedthroughs in the brass end 
plate. The complete assembly was inspected for dust and then placed into the sample cell. 

The sample cell had an outer diameter that was 38 mm, and an inner, cylindrical space 
was 38 mm long and 19 mm in diameter. One end of the cell was sealed by a sapphire 
window that had been coated with tin oxide to eliminate static charges. The other end w^s 
sealed by a brass plate containing five electrical feedthroughs connected to the oscillator 
and the four surrounding electrodes. The two sealing gaskets were made of indium-coated 
copper wire. The seals and feedthroughs were tested for leaks with a mass spectrometer 
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leak detector while the cell was pressurized with helium at 12 MPa. The cell was heated 
to 100 °C, pumped for four hours to remove volatile contaminants; then it was cooled to 
room temperature and isolated from the vacuum pump. The oscillator’s quality factor Q was 
monitored during the next four days. The time dependence of Q indicated that outgassing 
was negligible. 

The torsion oscillator’s resonance frequency / 0 , quality factor Q, and elastic aftereffect 
were measured in vacuum as functions of temperature. These measurements were used 
later in a model of the oscillator’s anelasticity to predict the frequency and temperature 
dependence of the oscillator’s spring constant and its internal losses [21] at all frequencies. 
Near T c , f Q ~ 11 Hz and Q ~ 1000. See Figure 6. No other mode was detected at frequencies 
between 0 Hz and 100 Hz. 

2.2. Loading xenon at critical density 

In order to obtain meaningful viscosity data close to T c , the density of the xenon near the 
oscillator had to be within a few percent of p c . To attain this we controlled both the average 
density and the density gradients. The average density was determined when the cell was 
loaded. 

After characterization of the oscillator in vacuum, the sample cell was filled with xenon 
(Matheson 99.995%) through a copper fill line. To adjust the loading, the cell was immersed 
and held horizontally in a thermostatted fish tank at a temperature less than 2 mK below T c . 
Xenon was added or removed to bring the liquid-vapor meniscus to the cell’s midplane. The 
fill line was crimped, cut, and sealed by soldering. Epoxy was poured around the electrical 
feedthroughs as a precaution against leaks. The cell was weighed one year before and one 
year after the flight. The difference between the weighings indicated that less than 0.02% of 
the sample was lost during this period. 

We observed the height of the meniscus when the sample was in equilibrium a few mK 
below T c . The inferred average density was (1 — 0.0015 ± 0.0017) p c . This difference from p c 
decreased the correlation length. At the lower bound of the fitting range, f = lx 10~ 6 , the 
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corresponding decrease of the viscosity at 0 Hz was (O.IOIq oo) %■ 

We measured T c by recording the temperatures at which the meniscus appeared and 
disappeared while the cell was immersed in the fish tank. The thermometer was a platinum 
resistance thermometer immersed near the cell. Measurements made three years and one 
year before the flight fell in the range T c = 289.736 ± 0.006 K. 

The xenon lowered the oscillator’s Q from 1000 to less than 1. This overdamping made the 
viscometer insensitive to vibrations associated with normal Shuttle operations and it led to 
other advantages. The overdamped oscillator was sensitive to viscosity changes in the range 2 
to 12 Hz. This allowed the oscillator to be calibrated by exploiting a hydrodynamic similarity. 
The data at multiple frequencies provided a powerful check on the viscometer’s accuracy, 
and they helped reveal the viscoelasticity of xenon near T c . Finally, the overdamping made 
the oscillator rugged enough to survive the strong accelerations associated with launch and 
touchdown. 

2.3. Thermostat 

Maintaining sample homogeneity near T c required that temperature differences within the 
sample be small. The CVX thermostat achieved such small temperature differences. It 
consisted of three concentric aluminum shells surrounding the thick-walled copper sample 
cell. The cylindrical shells and their end caps were made from 6-mm-thick aluminum with a 
radial gap of 13 mm between shells. The large radial gap and stiff, glass-filled polycarbonate 
spacers between shells made the design mechanically robust and insensitive to errors of de- 
sign and construction. The 38-mm separation between end caps allowed easy installation of 
the cell and its wiring. The weak coupling between shells that resulted from the large gaps 
increased the thermostat’s response time to more than one hour.' However, this was accept- 
able because the thermostat s response time was less than the sample’s internal response 
time near the critical point. The thermostat’s construction and operation are described in 
more detail elsewhere [18, 17, 22], 

The performance of the thermostat was verified by using a semiconductor thermopile to 
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measure the temperature difference imposed along a thin-walled steel sample cell both while 
controlling the thermostat at constant temperature and while ramping it at 10 /xKs -1 [22]. 
Extrapolating the results to CVX indicated that, while ramping at the slowest rate of -0.05 
/xKs -1 , the temperature difference along the thick-walled copper cell was only -0.11 /xK. 
We added half of this value to the difference caused by the thermistor’s power to obtain a 
maximum difference of 0.11 /xK of the cell’s wall temperature from its average value. This 
corresponded to a density difference of 0.003p c at t = 1 x 10~ 6 , which was acceptable. 

2.4. Flight cannisters 

The CVX flight package consisted of two sealed 0.8-m-taIl aluminum “Hitchhiker” cannisters. 
The “Experiment” cannister, shown in Figure 7, contained the thermostat and the more 
sensitive analog electronics. It also contained electrical batteries that were used to maintain 
the thermostat’s temperature above T c during descent when no power was available from 
the Shuttle. Keeping the sample above T c during descent prevented the formation of liquid 
whose sloshing could have damaged the oscillator. The “Avionics” cannister contained power 
conditioners, an accelerometer, and the data acquisition and control electronics. It also 
contained the four computers dedicated to the tasks of viscometry, temperature control, 
accelerometry, and communications. 

The cannisters were mounted in the Space Shuttle’s open payload bay, where the external 
heat load varied greatly depending on whether the bay was oriented toward deep space, the 
Earth, or the Sun. These heat load variations changed the Experiment cannister’s interior 
temperature, which in turn affected the gain of the viscometry electronics. To minimize 
these temperature changes, the sides of the cannister were thermally isolated and the mass 
of the cannister’s lid was increased to 45 kg. A heater attached to the underside of the lid 
provided additional temperature control. Between mission days 1 and 10, the Experiment 
cannister’s interior temperature was maintained within 1 K of 12 °C. 

The electronics and the thermostat were cooled only via radiation from the Experiment 
cannister’s lid. Thus, too large a heat load on the Experiment cannister would have caused 
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a disastrous loss of thermostat control. Covering the cannister lid with appropriate radi- 
ator tape and limiting the duration of the Shuttle’s Sun-facing orientations prevented this 
potential problem. 

2.5. Electronics 

2.5.1. Oscillator drive voltages 

The viscosity was deduced from the ratio of the torque applied to the oscillator to the 
deflection of the oscillator. Figure 8 is a schematic diagram of the circuit elements that 
were used to apply the torque and to measure the deflection. The paired 162 kfl resistors 
and 1000 pF capacitors allowed the simultaneous presence at the electrodes of the sub-audio 
frequency voltage used to drive the oscillator and the 10 kHz voltage used to detect the 
oscillator’s deflection. The 1 Mft resistor at the input of the lock-in amplifier grounded the 
oscillator at low frequencies. A smaller resistor would have decreased the signal-tonoise 
ratio unacceptably, and a larger resistor would have prevented effective grounding. The 
interelectrode capacitances C4, Cg, Cc, and C q shown in Figure 8 were approximately 0.3 
pF. This was much smaller than the cable capacitance C cable ] however, the 10 kHz signal 
was still detected with a satisfactory signal-to-noise ratio. 

In normal operation, the torque applied to the oscillator was proportional to a time- 
dependent voltage created by summing 400 equal-amplitude sine waves at frequencies evenly 
spaced from /, = 1/32 Hz to 12.5 Hz, 

400 

Vi,, (0 — A ^ sin [2nnf x t + <f> (n)] . (2.1) 

71 = I 

This waveform had a period of 32 s. The phases <fi (n) of the 400 components were chosen 
to minimize the waveform’s maximum excursion. A digital-to-analog converter with 16-bit 
voltage resolution and 512-Hz time resolution generated the waveform. A Iowpass filter 
smoothed the steps in the waveform. Occasionally, low-frequency measurements were made 
with a similar waveform whose duration was 512 s, and whose lowest frequency was /, = 
1/512 Hz. 
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Because the electrostatic torque was proportional to the square of the drive voltages, we 
used analog circuits to transform the input voltage Vi„. A bias voltage Vbc was added to 
Vj„ ( t ), and the square root of the sum, 

v, (i) = +V D0 [1 + (Kn (!) /Vdc)] 1/2 . (2.2) 

was obtained. This voltage was applied to one of the diagonally connected electrode pairs. 
A similar voltage, but of opposite phase, 

^2 (0 = -Vbc [1 - (M„ (t) /Vbc)] 1/2 , (2-3) 

was applied to the other electrode pair. (The opposite sign of the voltages V\ and V? reduced 
the torque’s nonlinearity.) The torque applied to the oscillator was approximately 

CaV dVosc 

Xq 

where L osc is the oscillator’s length, x ( t ) is the displacement at the oscillator’s tip, and 
Xq ~ 4 mm is the gap between the oscillator and one electrode. The second term of Eq.(2.4) 
represents a “softening” of the oscillator’s spring that lowered the oscillator’s resonant fre- 
quency / 0 . (See Section 4.3) The square-root circuit and the approximate symmetry of the 
electrode pairs made the torque on the oscillator very nearly a linear function of Vj n (£). The 
square-root circuit was necessary for linearity because the ratio V; n /Vbc was as large as 0.2. 
This ratio could not be decreased significantly by increasing the bias voltage above its actual 
value, Vbc = 30 V. Such an increase of Vbc would have risked pulling the oscillator against 
one of the electrode pairs, and it would have added nonlinearity to the oscillator’s equation 
of motion. The lowpass filter following the square-root circuit suppressed noise at 10 kHz. 

2.5.2. Oscillator displacement detection 

The oscillator’s displacement was detected by the unbalance of a capacitance bridge. See 
Figure 8. The bridge was driven by a 3 V rms , 10 kHz oscillator. An inductive voltage divider 
was adjusted to approximately balance the bridge. The out-of-balance signal was fed to a 


jUO , *(*) 

Vbc x o 


(2.4) 


13 


lock-in amplifier which generated a sub-audio frequency voltage K out that was linear in the 
difference of the capacitances 

AC ( t ) = (Ca + Cp) — ( Cs + Cq) . ( 2 . 5 ) 

This difference was approximately proportional to the oscillator’s displacement x(£). 

AC (t) ~ AC A • (2.6) 

The value of x$ was a compromise between increasing the sensitivity and increasing the 
nonlinearity as the gap was decreased. 

The viscometry’s input and output signals, V in and V(, ut were processed in five steps. (1) 
Anti-alias filters removed frequencies above 128 Hz. (2) The signals were simultaneously dig- 
itized at 512 Hz, in synchrony with the digital-to-analog converter which created V in ( t ). (3) 
The resulting time records were Fourier transformed, and all but the 400 lowest frequencies 
were discarded, thereby digitally filtering the data. (4) The transfer function, 

n ( t \ _ K»ut (/) 

^meas [J ) — -yT-Jjy, (2.7) 

was computed from the ratio of Fourier transforms of V in and V otA . ) (5) This function was 
stored as 401 complex numbers. In normal operation, the 32 s of data collection were followed 
by 32 s during which the data were processed, stored, and transmitted to ground. The same 
32 s waveform drove the oscillator during both halves of this 64 s cycle. The accuracy of the 
signal processing was verified by constructing a passive lowpass filter with a transfer function 
that resembled that of the overdamped oscillator. The filter’s transfer function measured by 
the CVX instrument agreed with that measured by a commercial spectrum analyzer (Hewlett 
Packard 35660A). It also agreed with the transfer function calculated from the values of the 
filter’s components. 

The inductive voltage divider shown on Figure 8 was developed by NIST’s Electricity 
Division to fit on a single circuit card. Such a small size was possible because CVX’s divider 
required only nine bits of resolution, much less than the twenty bits typical of commercial 
programmable dividers. The lock-in amplifier (Ithaco model 410) also fit on a single card. 
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2.5.3. Electric field effects 


Electric fields drove the oscillator. Thus, they were essential to CVX’s operation; however, 
they had two secondary effects. 

The first effect of the electric fields was an increase of the xenon density via electrostric- 
tion. Electrostriction was of greatest concern in the immediate vicinity of the oscillator be- 
cause the oscillator’s damping depended approximately upon a weighted integral of (r/p) 1 ^ 2 
over a volume within a viscous penetration length of the oscillator’s surface. Appendix B 
demonstrates that electrostriction had a negligible influence on CVX’s operation because it 
increased the average density near the oscillator by less than A p/p c = 0.001, the uncertainty 
in the sample’s average density. 

The second effect of the electric fields was electric-field-driven (dielectrophoretic) con- 
vection, which caused parcels of cooler, denser fluid to move toward regions of high electric 
field. Dielectrophoretic convection is analogous to the buoyancy-driven convection on Earth 
that transports cooler, denser regions that form near the top of a cell to the bottom of the 
same cell. In microgravity, such regions formed near the cell’s boundary when the cell was 
cooled, and they formed in the cell’s interior when the cell was warmed. Spinodal decom- 
position also caused their formation when the sample’s temperature was brought below T c . 
The characteristic pressure (chemical potential per unit mass) that drove dielectrophoretic 
convection in the gap between the oscillator and the electrodes was estimated a sPg~ 0.13 
mPa, four times greater than the hydrostatic pressure caused by microgravity. 

2.5.4. Oscillator amplitude effects 

During the experiment, the oscillator’s maximum displacement was 0.03 mm at the tip of 
the screen. The resulting product ■yr of shear rate 7 and fluctuation-decay time r was 
sufficiently low that CVX did not encounter near-critical shear thinning [23]. The oscillator 
dissipated approximately 7 pW in the xenon. The resulting rate of density change near 
the oscillator p was approximately proportional to the local power per unit volume Q 1 and 


15 



inversely proportional to xenon’s heat capacity at constant pressure Cp. The estimate 


PcCp 

integrated over the duration of the experiment, was negligible. 

Increasing the oscillator’s amplitude by a factor of two during ground 
that the viscometer’s response was independent of amplitude. 


( 2 . 8 ) 


tests demonstrated 


2.5.5. Temperature measurement and control 

Each thermostat shell had an embedded thermistor whose temperature was measured once 
every eight seconds with a resistance bridge. The shell’s temperature was controlled by a 
proportional-integral-derivative algorithm. The middle shell’s temperature was set 0.03 K 
below that of the inner shell, and the outer shell was set 0.3 K below that of the middle 
shell. A thermistor embedded in the cell’s copper wall operated at a power of approximately 
0.6 pW. Its temperature was read by an AC bridge and lock-in amplifier. 

By choosing each bridge’s reference resistor to have a value near that of its thermistor 
at T c , the need for an adjustable component, such as a ratio transformer, was eliminated. 
Temperature was inferred from the bridge’s unbalance instead of from the adjustment re- 
quired to balance the bridge. Far from T c , the gain of the lock-in amplifier was decreased 
to accommodate the bridge s large unbalance. Within 50 mK of T C) the rms scatter in the 
cell’s apparent temperature was approximately 10 pK. 

The cell’s thermistor was calibrated against the inner shell’s thermistor to simplify the 
cell s temperature control. The inner shell’s thermistor calibration consisted of a fit of the 
Steinhart-Hart equation [24] to the resistances at 0°C, 25°C, and 50°C. The manufacturer 
stated that these three calibration points had an uncertainty of 0.05 K. Thus, the uncertainty 
of the reduced temperature t was approximately 0.0024, and temperature scale nonlinearity 
contributed negligible error to the derived value of the viscosity exponent. As an independent 
verification of the calibration’s accuracy, the cell’s thermistor found T c = 289.721 K. To 
within 0.02 K, this value is consistent with the value T c = 289.736 ± 0.006 K measured 
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with a platinum resistance thermometer in the fish tank on the ground and with the value 
289.74 ± 0.02 quoted elsewhere for xenon [25]. The thermometry’s stability was verified by 
comparing the temperature of the cell T ce n with that of the inner shell T in . The difference 
Tceii — T; n measured during flight drifted less than 0.1 mK per day, and it differed from the 
difference measured nine months earlier by only 0.5 mK. 

3. Temperature timeline 

The density of a pure fluid near its critical point is extremely sensitive to temperature 
gradients. For example, at p c and T c + 1 mK, the isobaric thermal expansivity of xenon is 
more than a million times larger than that of an ideal gas; thus even a tiny temperature 
gradient can induce a significant density gradient. In the absence of gravity, this effect 
limits the fluid’s homogeneity. Once formed, the density gradient can be long-lived because 
microgravity allows very small values of thermal diffusivity Dt throughout the sample. At 
T c + 1 mK, the slowest thermal time constant calculated for the CVX sample in the absence 
of convection was about one week. 

CVX’s conservative design did not rely on convection to ensure that the density would 
be sufficiently uniform. Instead, temperature gradients were minimized by careful design of 
the sample cell, the surrounding thermostat, and the sequence of temperature changes, or 
“temperature timeline”. To keep the sample’s density acceptably close to the critical density 
Pc, the temperature timeline shown in Figure 3 used a two-part strategy. Far from T c , the 
temperature was changed by large, rapid steps, causing large, temporary inhomogeneities in 
the sample. Each step was followed by a waiting period which exceeded the xenon’s longest 
equilibration time constant and which brought the xenon’s density close to p c . Close to T c , 
although the xenon’s temperature was changed without waiting for equilibrium, the density 
remained sufficiently close to p c to obtain meaningful measurements of viscosity. 

Candidate timelines were tested with a numerical model of entropy transport within the 
sample. The model contained two simplifications which allowed efficient testing. First, the 
xenon sample was modeled as an infinitely long cylinder whose density depended only on t he 
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radial coordinate. In this approximation, the model was one-dimensional, and heat conduc- 
tion through the cell’s internal parts was ignored. Second, xenon’s properties were estimated 
by an approximation to the cubic model equation of state. (The model is summarized in 
Reference [2].) In this approximation, terms of order 6 2 and higher were dropped to remove 
the need for iterative calculations. Keeping the state parameter 0 < 1 made the approxi- 
mation valid. Figure 9 shows the density deviations calculated for a timeline similar to that 
used by CVX. Below T c + 100 mK, the density deviation in the cell’s interior (r /Rc eil = 0) 
is less than 0.13%. Heat conduction through the copper wires that supported the electrodes 
reduced this to less than 0.06%. 

Our confidence in the model s physics came from a recent microgravity experiment [4], 
in which Wilkinson et al. demonstrated agreement between the measured and calculated 
values of thermal equilibration time constants in a sample of SFg near the critical point. 

Our confidence in the numerical calculations came from a numerical calculation that 
could also be described by two analytic calculations. The numerical example calculated the 
density change A p {t) in the sample’s interior during and after a 2000 s ramp of the boundary 
temperature from T c + 0.8 K to T c + 1.0 K. During the ramp, A p {t) increased until it reached 
its maximum value at 2000 s. After the ramp, A p (t) decayed, becoming exponential after 
8000 s. The decay s time constant agreed with the first analytic calculation, which gave 
5200 s for the cylinder’s slowest radial mode. The maximum value of 0.0022p c agreed with 
the second analytic calculation. This calculation was the adiabatic upper bound to the 
density change which occurs in the cell’s interior following a sudden change of the cell wall’s 


temperature from T, to T 2 . The bound is 

rTj 




WPJ 


r'J 2 

J 


c v d.T, 


(3.1) 


(TJP C ) ( dP/dT ) pc J Ti 

where c v is the constant- volume heat capacity. Eq.(3.1) shows that changing the temperature 
from T c + 0.1 K to T c induced changes of density in the sample’s interior that were less 


than 0.2%. Because this constraint is independent of the rate of temperature change, the 
maximum temperature ramp near T c was limited by other considerations. The ramp rate 
had to be slow enough that the viscosity would not change significantly during the time 


18 



required to measure it. 

The timeline used by CVX included the following features. 

1. Each temperature step was followed by a waiting period at constant temperature to 
observe the sample’s approach to equilibrium. 

2. The viscometer’s calibration data were taken during the waiting period at T c + 1 K. 

3. An initial “fast” temperature ramp passed through T c at the rate -1 /i K s _l and located 
T c on the cell’s thermometer to within 0.1 mK. 

4. A later series of “slow” temperature ramps, the slowest of which passed through T c at 
the rate -0.05 /iK-s" 1 and collected most of the data near T c . 

4. Data reduction 

During the mission, most data were downlinked in nearly real time for preliminary analysis, 
thereby allowing adjustments to the timeline. After the mission, all of the experimental data 
were retrieved from CVX’s hard disk. They comprised approximately 10 4 measurements of 
the transfer function G meiis (/), plus accompanying temperature measurements. Figure 10 
shows typical measurements of magnitude and phase of G meas (/). After culling and averaging 
the data, we used a portion of the data from T c + 1 K to calibrate the viscometer. Viscosity 
data were then derived from the averaged G meas (f) by combining the calibration with the 
model for the oscillator’s equation of motion. 

After the mission, data farther above T c were collected. For these ground measurements, 
we calibrated the viscometer at T c + 5 K instead of at T c + 1 K due to concerns about internal 
waves at low frequencies near T c [26]. Appendix A lists viscosity values derived from the 
microgravity and ground measurements. 
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4.1. Culling and averaging of the data 

The orbital environment influenced CVX in unexpected ways. Figure 11 shows |G mea5 (1 Hz)| 
measured during a 24-hour period. The data include oscillations with a 45-minute period and 
three spikes which lasted approximately five minutes each. Such oscillations and spikes 
occurred throughout the mission, and they affected the magnitude of the transfer function 
at all frequencies. The oscillations period was half of the Shuttle’s orbital period, which 
suggests that they were not driven directly either by the Shuttle’s exposure to the sun or by 
the direction of the Earth’s magnetic field. In ground tests, exposure of CVX to magnetic 
field variations comparable to those in orbit did not influence G meas (/). The minima of the 
oscillations occurred when the Shuttle was near the equator, and the spikes usually occurred 
when the Shuttle was near the South Atlantic Anomaly, a region of minimum magnetic field 
near Argentina. For low Earth orbit, both regions are associated with larger fluxes of the 
charged particles trapped in the Earth’s magnetic field. Thus, the timing of the oscillations 
and the spikes, and their absence in ground tests of CVX, are consistent with the hypothesis 
that they were caused by charged particles. 

Because the oscillations and the spikes were present only in the magnitude data, they were 
a time-dependent response of the instrument and not of either the sample or the oscillator. 
We excluded the spikes from the data, and we suppressed the oscillations by averaging 
the remaining data in groups of 45 minutes. Each group was thus the average of up to 42 
measurements. We did not weight the averaged data by the number of measurements in each 
group. Instead, we simply excluded from groups containing fewer than 30 measurements. 
The cell s temperature was averaged in a corresponding manner. 

Figure 12 shows that the orbital environment also greatly increased the noise at frequen- 
cies below 0.1 Hz. In contrast to the oscillations and the spikes, this very low-frequency 
noise was present in both the phase and the magnitude, suggesting that either the sample, 
the oscillator, or both were directly affected. This in-orbit noise prevented our use of the 
low frequency measurements for deriving the transducer factor k lr /k 0 (see Section 4.3). 

In addition to data adversely affected by the orbital environment, we excluded data taken 
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while the sample’s density distribution was far from equilibrium. Such data were taken below 
T c and during and following large steps of temperature. 


4.2. Hydrodynamic similarity and the viscometer’s calibration 


The viscometer’s calibration exploited a hydrodynamic similarity which applies to the hy- 
drodynamics of an immersed oscillator [27]. This analysis assumed that the fluid was homo- 
geneous and that the oscillator moved as a rigid body. The similarity can be introduced by 
considering the transfer function of a damped harmonic torsional oscillator. The oscillator’s 
response to a steady sinusoidal force at frequency / is given by 



(4.1) 


where k, Q , and / 0 are the oscillator’s spring constant, quality factor, and undamped res- 
onance frequency respectively, and G (/) is the ratio of the angular displacement of the 
oscillator to the external torque applied to the oscillator. By giving Q a frequency depen- 
dence appropriate for linearized hydrodynamics, this expression can describe the response 
of an immersed oscillator over a wide range of frequencies. For small amplitude oscillations, 
the response is 

'■ <42> 

where R and p s are a characteristic length and density of the oscillator respectively. The 
function B ( R/S ) characterizes the oscillator’s damping in terms of the viscous penetration 
length, defined by 



where rj and p are the fluid’s viscosity and density. At the nominal frequency of 5 Hz, 
6 ~ 60 pm. Thus 6 is larger than the nominal diameter of a screen wire (2 R = 28 pm), but 
much less than the distance between wires (847 pm). 

The viscosity enters G (/) only through the function B ( R/6 ). Most oscillating-body 
viscometers are sensitive to viscosity changes only in a narrow range of frequencies near / 0 , 
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and B{R/6 ) must be calculated from the hydrodynamic theory for an idealized geometry. 
In contrast, we determined B ( R/6 ) by measuring G (/) when the xenon was in a reference 
state of known tj and p far from T c and then inverting Eq.(4.6) below. This procedure 
calibrated the viscometer for changes in the viscosity with respect to the reference viscosity. 
We emphasize that this calibration did not require knowledge of the geometry of either the 
oscillator or the surrounding electrodes, and it was indifferent to the choices of the parameters 
R and p„. The accuracy of this technique was demonstrated in Reference [27]. There, the 
viscosity of carbon dioxide was measured using a viscometer much like the CVX viscometer, 
and the results agreed with accurate published data acquired by conventional means. 

The CVX viscometer’s calibration was derived during flight from an average of four 
45-minute measurements of G (/) near the temperature T c + 1 K. At this temperature, we 
equated the viscosity with that determined by the high-Q torsion oscillator [10], which had an 
uncertainty of 0.8%. After applying the corrections described below, useful values of B {R/6) 
were derived from the transfer function at frequencies from 0.3 to 12.5 Hz, corresponding to 
a range of 6 much broader than that encountered upon approaching T c at fixed frequency. 
We represented the magnitude and phase of B {R/6) by polynomial functions of In {R/6) in 
the range 0.05 < R/6 < 0.4. 

Figure 13 shows the measured values of B {R/6) and an approximation to B {R/6) derived 
from an analytical model of the oscillator. This model is Stokes’ solution for the viscous force 
exerted on a transversely oscillating cylinder. We used the modern formulation [28], 


B C y 1 {R/6) = -i 


, 4 AT, (z) 1 
zK 0 (z ) 1 


where K n {z) is the modified Bessel function of order n, and z 


(l + i){R/6). 


(4.4) 


4.3. Corrections to the ideal transfer function 

Figure 10 shows an example of the measured transfer function C meas (/). Several effects 
caused G meas (/) to differ from the response G (/) given by Eq.(4.2). The largest effect was 
due to the driving and detection electronics shown in Figure 8, which acted like a lowpass 
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filter with a transfer function G^. The transfer function of the electronics was obtained 
by disconnecting the driving and detection stages from the cell and measuring their transfer 
functions independently. 

The oscillator’s torsion spring was not ideal because slow relaxation processes within the 
nickel itself caused internal friction and creep. At the low amplitudes used in our measure- 
ments, these phenomena could be described as linear anelastic effects [21]. We corrected for 
anelasticity by generalizing the spring constant A; to a complex function of frequency 


k -Ml + (/)+#(/)], 


(4.5) 


where k 0 is the real part of the spring constant at the resonance frequency. The real functions 
ijj (/) and 4> (/), were measured in vacuum from 0.001 to 10 Hz at several temperatures [21]. 

The oscillator’s resonance frequency / 0 was corrected for its dependence on temperature 
and for the electrostatic spring softening induced by the DC voltage Vdc applied to the 
cell’s electrodes. Both effects could be measured in vacuum only. Due to xenon’s dielectric 
constant, we expected a 0.4% decrease of / 0 upon filling the cell, and this was included in the 
correction to / 0 . (The function B(R/6 ) accounted for the xenon’s effect on the oscillator’s 
hydrodynamic mass.) 

We did not correct G meas (f) to account for a second mode of the oscillator at higher 
frequency. A mode near 54 Hz, near the calculated bending frequency, was seen in vacuum 
measurements of a similar oscillator. However, no such mode was seen below 100 Hz with the 
CVX oscillator. Perhaps the bending mode was shifted to higher frequencies by stiffening 
caused by a slight crease of the nickel screen. 

With these corrections, the transfer function becomes 


G m eas{f) ~ G e u* (/) 


Wr 


1 4 V (/) + i* (/) - (jJ + t f£ J M B (R/6) 


-1 


where k tr is a real constant given by the product 

/ oscillator torque \ / output voltage \ / 1 
>r \ input voltage J \ oscillator displacement/ \G e | 


(4-0) 


(4.7) 
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The product of transducer factors k tr related the measured voltages V in and V out to the 
oscillator’s torque and displacement. We assumed that k tr was independent of frequency. 

Determination of the ratio k tr /k 0 requires accurate knowledge neither of k 0 nor of the 
transducer factors. For example, in ground tests we obtained kt r fk^ from measurements of 
Gmeasif) made near 0.01 Hz, where viscous damping was insignificant. Derivation of the 
value 

ktr ICn^O-Ol Hz) | 

k 0 1- -0(0.01 Hz) (4 - 8 ) 

required knowledge only of the real part of the anelastic correction. In orbit, we were unable 
to measure k tr /k 0 in a similar fashion because the signal-to-noise ratio of G meas (0.01 Hz) 
decreased by a factor of 10 or more. Instead, we set k tr /k 0 to the value that made Re (77) 
independent of frequency for / > 1 Hz in the temperature range 3 x 10 -5 < t < 5 x 10 -2 
where viscoelastic effects were absent. The in-orbit and ground values of kt r /ko agreed within 
0.3%. 

The parameter kt r /kQ had linear dependences on the Experimental cannister’s temper- 
ature, the cell’s temperature, and time. The effects of these dependences can be seen in 
Figure 14, which is a plot of | G meas (0.5 Hz)|. At 0.5 Hz, \G meas \ is nearly independent of the 
sample s viscosity; thus, only unintentional dependences are present. The dependence on the 
Experimental cannister’s temperature probably originated in the electronics; its coefficient 
was determined in ground tests. The dependence on the cell’s temperature was assumed to 
result from the temperature dependence of the elastic constant of the torsion fiber and the 
thermal expansion of the capacitors. The torsion constant’s temperature dependence was 
inferred from measurements of df qJ dT } and changes in the capacitance were determined by 
fitting |G meas (0.5 Hz) | to a linear function of the cell’s temperature. The coefficient of the 
linear time dependence was determined by fitting to |G meas (0.5 Hz) |. The origin of this time 
dependence was unknown. However, a similar time dependence occurred in |C7 meas (0 Hz)|. 
which was linear in the oscillator’s average displacement; both time dependences may have 
resulted from creep in the oscillator’s equilibrium position. 
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(T - Teal) raised to a small power such as z v , all of our results would have been unchanged. 

As an additional check on the data reduction, we systematically varied various data 
reduction parameters from their nominal values. We changed the density by 1% from p c . 
We changed the calibration viscosity by 3%. We changed the oscillator’s vacuum resonance 
frequency by 3% from / 0 . We perturbed the transfer function G e iec used to describe the 
electronics. None of these variations caused significant changes in the data’s consistency. 


5. Results 

5.1. Sample homogeneity 

5.1.1. History independence 

Density inhomogeneities that were induced by temperature changes caused no problems. 
Figure 15 supports this claim by showing the oscillator’s phase measured at 2 Hz during a ten- 
day period. On Figure 15, measurements made at 300 mK, 100 mK, 30 mK, and 3 mK above 
T c are highlighted by boxes. At each temperature, measurements were repeated as much as 
seven days later and with very different temperature histories. The phase measurements 
agreed to within 1 mrad. This corresponds to an agreement in density of 0.1%. 

Figure 4 shows that the sample’s density could be kept near p c even within 0.1 mK of T c . 
Data collected during three ramps through T c are shown. The slowest ramp, at -0.05 pKs 1 , 
collected most of the data near T c . The medium ramp, at -1 pK s 1 , was intended only to 
locate T c to within 1 mK. Nevertheless, it gave results above T c that were indistinguishable 
from those of the slowest ramp. Furthermore, after culling data points as desciibed eailiei, 
the data at higher temperatures also were independent of the sample’s temperature history. 

On Figure 4, data from the fastest ramp at -33 pKs~ l are displaced above and to the 
left of the equilibrium data. Because this ramp started at T c + 0.2 K, Eq.(3.1) permits 
the sample’s interior density to decrease by up to 0.3%. This decrease is consistent with 
the data’s vertical displacement. (Also, the thermometer’s lag of several seconds caused a 
horizontal displacement of 0.1 mK.) Note also the ramp-rate dependence of the phase at 2 
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Hz below T c . This time dependence is consistent with phase separation. 

5.1.2. Equilibration 

When the cell’s temperature was held constant, the viscosity rapidly approached a steady 
value, typically within one hour. (See Figure 15.) Even after accounting for heat conduction 
through the cell’s metal electrodes, this was approach was faster than expected from the 
xenon’s thermal diffusivity. At 300 mK, 100 mK, 30 mK, and 3 mK above T c , the sample’s 
calculated time constants were 1 h, 2 h, 6 h, and 30 h, respectively. 

Three phenomena contributed to the rapid approaches to equilibrium. First, convection 
driven by gravity and by electric fields slowly stirred the xenon. The stirring was effective 
because the open geometry around the oscillator reduced viscous damping of the convective 
flow. For example, if a density deviation A p = 0.01p c had occurred in a region of charac- 
teristic size L — A mm, an acceleration 10~ 6 of Earth’s gravity would have turned over that 
portion of the sample in a time of order ? ?/ (gLAp) « 3 minutes. Second, the oscillator was 
thermally isolated from the cell wall, and it had negligible heat capacity. Thus, the deviation 
of the density caused by a temperature step above T c was equal to that of the cell’s interior 
(Eq.(3.1)) and not of the much larger deviation that occurred near the cell’s boundaries. For 
example, the initial density deviation in the cell’s interior following the steps from T c 4- 300 
mK to T c + 100 mK was only 0.3%. Third, the oscillator was insensitive to deviations of 
the density from p c if the deviations averaged to zero over the oscillator’s surface. This 
insensitivity occurred because the oscillator’s response was approximately proportional to 
the surface-averaged value of Jgp and, above T c 4- 1 mK and for small deviations from p c , 
the viscosity is proportional to the density. 

To demonstrate the viscometer’s sensitivity to density inhomogeneities, we deliberately 
created an inhomogeneous sample by cooling the cell to T c - 9 mK and then warming it to 
T c +0.5 mK. We then observed the apparent viscosity for four hours. This set of observations 
is labeled ‘B” and enclosed by an oval on Figure 15. It can be compared to the set labeled 
A that was cooled to T c - 0.5 mK and then warmed to T c 4-3 mK. Because xenon’s thermal 
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diffusivity at T c + 0.5 mK is three times smaller than at T c + 3 mK, one might expect the 
approach to equilibrium of set B to be three times slower than the approach of set A. Instead, 
the approach of set B was slower by much more than a factor of three; the apparent viscosity 
differed from the equilibrium viscosity by 5% at the end of the four-hour observation. 

One contribution to set B’s slower equilibration was the limited thermal conductance 
of the copper wire that connected each electrode to the cell wall. Below approximately 
T c + 1 mK, the conductance of the copper wire was exceeded by that of the xenon between 
the electrode and the cell wall, and the length //thermal that governed the sample’s thermal 
diffusion was the cell’s 19 mm diameter. In contrast, //thermal for set A was the 8 mm gap 
between the electrodes. 

Another possible contribution to set B’s slower equilibration was the size of its inhomo- 
geneities. They were formed at T c — 9 mK, where surface-tension-driven flows ensured that 
the liquid and vapor regions quickly grew as large as the sample’s dimensions. In contrast, 
set A’s inhomogeneities were formed between T c and T c — 0.5 mK during a two-hour period; 
this is sufficiently brief that the liquid and vapor regions may not have had enough time to 
grow as large as the 8 mm interelectrode gap [7]. In that case, //thermal was limited to the size 
of set A’s inhomogeneities, thereby speeding the decay of the inhomogeneities by thermal 
diffusion. 

5.2. Fits to the data 
5.2.1. Model of the data 

Our model of the frequency-dependent data was Eq.(1.2). The frequency scaling function 
S ( 2 ), including corrections to two typographic errors in Ref. [8], is given in Appendix C. We 
introduced the parameter A into the argument of the scaling function to obtain agreement 
with our data. Setting A = 1 caused large systematic deviations from the fits. 

In Eq.(1.2), the background viscosity was described by the sum of analytic functions of 
the density and the temperature: 

Vo(T,p) = Vw{T) +Vo\ {p)- (5-1) 
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These functions were determined from the viscosity measured by others far from the critical 
point. (See Appendix C of Reference [10].) For consistency with the measurements of 
Reference [10], we set -q 0 (T c ,p c ) = (51.3 ± 0.4) /iPas. The present results are not sensitive 
to this description of 770. We used the value p c = 1116 kgm -3 determined from Narger and 
Balzar ini’s Sample No. 1 [29]. Both Sample No.l and CVX used high purity xenon from 
Matheson. We estimated the uncertainty of p c as the difference of 1.3 kg m -3 between the 
values for Samples No. 1 and 2. 

Ideally the exponent characterizing the divergence of any thermophysical property would 
be determined from data so close to T c that the critical anomaly is much larger than the non- 
critical background. The small value of the exponent z n makes this impossible for viscosity, 
and the most accurate determination of z rt from the data requires knowledge of the back- 
ground viscosity 770 and of the crossover function H (£). Recently, Luettmer-Strathmann, 
Sengers, and Olchowy [30, 31] published equations designed to describe the crossover behav- 
ior of both viscosity and thermal conductivity with a single fluid-dependent wave-number 
cutoff q D . We used the older crossover function by Bhattachaxjee et al. [9] for viscosity 
because it is simpler and has the same asymptotic behavior as the newer function. The older 
crossover function makes approximations not used in the newer function. The chief approx- 
imation is neglect of the divergence in the constant- volume heat capacity c v [32]. Because 
we fitted to data only at reduced temperatures less than 10~ 3 , where the heat capacity ratio 
c v /cp < 0.002, this approximation introduced negligible errors. 

The following subsections describe the results for T c , z v , A , q c , and q D obtained by fitting 
the model to the viscosity data. The fit included the values of Re (77) and Im (//) measured 
at 2 Hz, 3 Hz, 5 Hz, 8 Hz, and 12 Hz within the range 10~ G < t < 10 -4 . Expanding the 
range of fitted data to 10~ 6 < t < 10“ 3 caused negligible changes in the fitted parameters, 
while slightly increasing the fit’s systematic deviations: The value of reduced chi-squared 
increased from 0.97 to 1.05. 



5.2.2. Viscosity cusps near T c 

The fit ignored the data within 0.3 mK of T c and yielded T c = (16.571237 ± 0.000008) °C. 
The standard uncertainties quoted here and below represent a coverage factor of one (68% 
confidence interval), and they allow for the correlations among the parameters. (We recall 
that our thermometer’s uncertainty with respect to ITS90 was 50 mK.) The phase and 
magnitude data exhibited cusps close to the fitted value of T c . The minimum in the phase 
at 2 Hz shown in Figure 4 is an example of such a cusp. As illustrated in Figure 16, the 
phase at frequencies higher than 2 Hz exhibited a sharp local maximum. Table 1 lists the 
temperatures of three example cusps obtained by separate fits to data in a 2 mK span that 
included T c . 


Table 1: Cusps near T c 

Cusp 

(T cusp - T c ) / fiK 

phase at 2 Hz 

+56 

magnitude at 5 Hz 

-1 

phase at 5 Hz 5 Hz 

-17 


Each value of T cusp was determined to within 10 /xK. The values of T cusp are thus mutually 
inconsistent even though they all differ from T c by less than 60 fiK. We allowed T c to be a 
fitted parameter because T c ^ T cusp . Viscoelasticity and time-dependent phase separation 
determine T cusp . 

5.2.3. Viscosity exponent z v 

The viscosity exponent was obtained by dividing the fitted value, = 0.04349 ± 0.00035, 
by the correlation length exponent, v — 0.630, [33, 34] to obtain z n = 0.0690 ± 0.0006. This 
value for z n is 7% larger than the value obtained from the high-Q torsion oscillator [10], even 
though the two data sets agree at t > 3 x 10 — 4 , where gravity had negligible influence on the 
high-Q data, and even though both data sets were analyzed with the same crossover function. 
The high-Q data led to a lower value of z n because they were obtained at larger reduced 
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temperatures, where the inaccuracy of the crossover function used here and in Ref. [10] is 
more important. Consistent with this attribution, fits to the present data at larger reduced 
temperatures yielded smaller values of z v . These fits also showed systematic deviations that 
would not have been detected in the high-Q data because the latter had more noise. 

Two recent papers by Hao, Ferrell, and Bhattacharjee have refined the theoretical esti- 
mate of Zq. The first paper [35] is a mode coupling calculation of z v in a two-term epsilon 
expansion. In contrast to earlier work by Siggia et al. [36], Hao et al. found that vertex 
corrections were significant, and they argued that a two-loop perturbation expansion in three 
dimensions would give a more accurate result. The second paper [15] gives the more accurate 
calculation, which leads to z n = 0.067 ± 0.002, where the uncertainty is the authors’ “rough 
estimate of the inaccuracy introduced by various approximations”. Thus the experimental 
and theoretical values of z v agree within their combined uncertainties. 

5.2.4. Viscoelastic time scale A 

The value of A was obtained by replacing the fluctuation-decay time r 0 with the free pa- 
rameter Tq. This procedure yielded Tq = (2.31 ±0.06) ps, thereby determining the fre- 
quency scale Az with a relative uncertainty of 0.03. However, the relative uncertainty of 
the ratio A = Tq/t 0 = 2.0 ± 0.3 is 0.15 because of the larger uncertainty of the fluctuation- 
decay time amplitude r 0 , given by Eq.(l.l). For xenon, r 0 = (1.15 ±0.17) ps, and the 
uncertainty of r 0 is dominated by the uncertainty of xenon’s correlation length amplitude, 
£o = (0.184 ± 0.009) nm, which we estimated as the difference between the two measurements 
reported by Giittinger and Cannell [37], 

Bhattacharjee and Ferrell’s theory of near-critical viscoelasticity [38, 8] was stimulated 
by Bruschi et al .' s high-frequency measurements of viscosity near the critical point of CO 2 
[11]. These measurements, and later measurements by Bruschi in CO 2 [12] and by Izumi 
et al. in xenon [13], showed only qualitative agreement with theory. Berg and Moldover’s 
single-frequency measurements of two binary liquids required viscoelasticity in the analysis. 
(See Figure 9 in [14].) However, the magnitude of the observable viscoelastic effect was only 
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0.5%, in hindsight we see that their data are consistent with A > 1 . 

In contrast to these previous measurements, the present data provide a quantitative test 
of the predicted functional forms. As shown by Figure 2, the functional form S ( Az ) describes 
the data for both Re ( 77 ) and Im ( 77 ) over a wide range of scaled frequency (0.0004 < Az < 33), 
provided that the viscoelastic time scale is adjusted by the factor A = 2.0±0.3. In particular, 
the description is consistent with the prediction that Im ( 77 ) / Re ( 77 ) approaches 0.035 for large 
values of z. 

5.2.5. Wavevectors q c and q D 

Fitting to the logarithms In {q c £ 0 ) and In (q D £ 0 ) instead of q c £ 0 and q D £ 0 made the fitting 
routine more robust. The values obtained were qdo = 0.051 ± 0.007 and q D ^ 0 = 0.16 ± 0.05. 
The fitted value of q c (o agrees with the independent value 0.059 ± 0.004 determined from 
published data that do not include the near-critical viscosity. See Appendix D. 

6. Acknowledgements 

We thank the CVX development team led by A.M. Peddie and J.L. Myers at NYMA Inc., 
the TAS -01 Mission team led by N.F. Barthelme, and the STS-85 crew. We thank R.A. 
Ferrell, R.W. Gammon, and R.A. Wilkinson for stimulating and critical comments. The 
CVX project was managed by R.W. Lauver and I. Bibyk through the NASA Lewis Research 
Center and was funded by the NASA Microgravity Science and Applications Division. 

A. Appendix: Tablulated viscosity data 

Table Al and A 2 give values of Re ( 77 ) outside the viscoelastic region. Only data taken at 2 
Hz are listed. The data at higher frequencies were consistent with those at 2 Hz; however, 
they had smaller signal-to-noise ratios. Some of the values represent a single 45 -minute 
data point, and their uncertainty was estimated as the standard deviation of all of the 
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45-minute points from a fit to the data in the range 10~ 6 < t < 10 -3 . The other values are 
averages of N 45-minute points, and their uncertainty was estimated as ax/y/N. 

Table A3 gives values of Re (77) and Im (77) / Re (77) inside the viscoelastic region. Each 
temperature represents a single 45-minute data point. 
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B. Appendix: Electrostriction 


The relative increase in the density A pj p c caused by a time-independent, electric field E is 


A £ 

Pc 



(B.l) 


where Pe — (cxeAo — 1) £qE 2 /2 is the induced change in chemical potential per unit mass, 
P c is xenon’s critical pressure, and cq and ex e are the dielectric constants of vacuum and 
xenon. Electrostriction is conspicuous near T c because of the large density changes resulting 
from the large value of Xt> the reduced susceptibility at constant temperature. 

As a first step toward calculating E , we ignored the oscillator’s presence and modeled the 
electric field between the viscometer’s driving electrodes as that in a parallel plate capacitor. 
See Figure 17. If the voltages on the two electrodes are V\ and V 2} and the distances from 
the oscillator to the two electrodes are X\ and x 2 , this electric field is 


(Vi ~ Vi) 

(Xj + x 2 )' 


(B.2) 


Electrostriction due to E 0 was insignificant. 

Next we modeled the field concentrated near one of the thin wires comprised by the 
oscillator. We adapted Spangenberg’s [39] model for the electric potential in a vacuum tube, 
consisting of two electrodes and a grid modeled as planar array of parallel line charges. The 
radial field E r near one of the wires, modeled as a cylinder, is 


E r (r,6) = (E l +E 2 


/ a \ ^ 

R \ 2 

( h — ) Eo 

1 + - 

\47t r) 

W 


cos (6) . 


(B.3) 


Here r and 0 are cylindrical coordinates centered on the wire, R is the radius of the wire, 
and a R is the spacing between wires. The first term in Eq.(B.3) originates in the wire’s 
induced line charge, and it disappears when the fields E\ = V\/x\ and E 2 = V 2 /x 2 are equal 
and opposite. It is a factor of two smaller than in Spangenberg’s model because the charge 
induced on the CVX oscillator was spread over transverse as well as longitudinal wires. The 
second term originates in the wire’s induced line dipole; it was added to Spangenberg’s model 
to create a zero potential surface at r = R. 
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We applied Eq.(B.3) to CVX using the values a = 0.85 mm, R = 0.009 mm, and 
V\ = ~V 2 = V DC = 30 V. We allowed for asymmetry in the viscometer’s construction by 
estimating the unequal distances i,=3mmandi 2 =5 mm. At the wire’s surface, the field 
was concentrated by the factor E r /E 0 = 6. At T c + 0.3 mK, the associated density increase 
was only 0.05% and thus negligible. 

C. Appendix: Frequency-dependent scaling function S ( z ) 

Bhattacharjee and Ferrell calculated the frequency-dependent viscosity of a classical fluid 
near its critical point. The relevant portion of their results in Reference [8] are reproduced 

here, correcting two typographic errors in the original published expressions for ln(5' 4 ) and 
for R (z). 

The dependence of the viscosity on correlation length £ and frequency f is 

v (£, f) = v (f, o) s z ” / ( 3+2 ’d, (c.i) 

where the argument of the scaling function S (z) is the scaled frequency defined by 2 = 
-mfr. Bhattacharjee and Ferrell used the decoupled-mode theory to calculate S (z) to 
single-loop order. They accurately approximated their result in closed form by an average 
of calculations in two and four dimensions. 

S^(s 2 y'\s,f , 02 ) 

Here, the tilde refers to additional rescalings given by 

S 7 {z) = S 2 ((2/e) 3 z) and 5 4 = S 4 ((8/e) z ) . (C.3) 

The scaling functions S 2 and S 4 are given by 

M *)- (rrT^.GrrT)’ 1 "* (^) 

and 

MS 4 )= (“ 3 + f) + ( 1 ~ f + p?) ln - + + (C o) 
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E. Figures 


1. Log-log plot of xenon’s viscosity measured near the critical point. The asymptotic line 
has the slope z^/u = 0.0435 deduced from the present microgravity data. Near T c , the 
CVX microgravity data (Re (77) at 2 Hz) depart from the asymptotic line because of 
viscoelasticity. The two sets of ground data depart from the asymptotic line further 
from T c because the xenon stratified in Earth’s gravity. 

2. Xenon’s viscosity at critical density measured at frequencies from 2 to 12 Hz. The solid 
curves resulted from fitting Eq.(1.2) to the data in the range 10 -6 < t < 10 -4 . (a) The 
real viscosity Re (77). Near t = 10 5 , the data depart from the 0 Hz curve because of 
viscoelasticity, (b) The ratio Im(77 ) /Re ( 77). For clarity, the ratio data at frequencies 
above 2 Hz are displaced downward by integer multiples of 0.005; otherwise they would 
coincide at t > 10 -5 . 

3. The temperature timeline followed by CVX. Both panels show the same data but with 
different vertical scales. 

4. Viscometer response as a function of temperature very near T c . Both panels show data 
at three ramp rates. On the lower panel, the temperature scale is expanded by a factor 
of 8. (a) Above T c , data collected at the fastest ramp rate of -33 /zK-s -1 reflect the 
0.3% lower density caused by starting the ramp at T c + 0.2 K. (b) Above T c , data 
collected while ramping at -0.05 /zKs _1 agree with those collected five days earlier 
while ramping 20 times faster. 

5. Cutaway view of the CVX viscometer cell. The cylindrical volume occupied by the 
xenon was 38 mm long and 19 mm in diameter. Torque was applied to the screen by 
applying different voltages to diagonal pairs of electrodes while maintaining the screen 
at ground potential. 
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6. The magnitude of the oscillator’s transfer function measured in vacuum and in xenon 
at critical density. 

7. Cross-section of the Experimental cannister. 

8. Simplified diagram of the circuit used to drive and detect the oscillator’s motion. 

9. The normalized density deviations calculated for a timeline similar to that used by 
CVX. Each curve was calculated for at a radius r within a cylinder with isothermal 
walls at radius Rc e \\- 

10. The oscillator’s transfer function G meas (/) in the frequency region used to determine 
the viscosity. 

11. Oscillations and spikes in the magnitude data at 1 Hz. The oscillation’s 45-minute 
period was half that of the Shuttle’s orbital period. The minima of the oscillations 
occurred when the Shuttle was near the Earth’s equator. The three spikes indicated 
by arrows occurred when the Shuttle was near the South Atlantic Anomaly. 

12. Examples of the transfer function’s magnitude measured at very low frequencies, where 
noise was much greater in orbit than on Earth. 

13. The calibration function B(R/8). The smooth curves indicate B(R/6) calculated 
from the hydrodynamic theory of a transversely oscillating circular cylinder of radius 
R = 13.4 /im. (a) Magnitude scaled to reveal departures from the dominant (R/6) 3 ^ 2 
behavior, (b) Phase. 

14. The time dependence of the oscillator parameter k tr . The solid curve is a fit to the 
transfer function’s magnitude at 0.5 Hz, assuming that k lr was a linear function of the 
cell’s temperature, the Experiment cannister’s temperature, and time. 

15. The transfer function’s phase at 2 Hz as a function of time. The height of the boxes 
corresponds to 0.2% of the density. The reproducibility of the oscillator’s phase at 300 
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mK, 100 mK, 30 mK, and 3 mK above T c demonstrates that density deviations could 
be reduced to within 0.1% of the average density at these temperatures. Filled points 
were obtained below T c . The points denoted “A” and “B” were obtained shortly after 
warming the sample from below T c . As discussed in the text, the equilibration of set 
B was much slower than that of set A. 

16. Cusps in the data very close to T c . These data were taken at 5 Hz; similar features 
occurred at other frequencies above 2 Hz. (a) Maximum damping (minimum ampli- 
tude) of the oscillator occurred within 0.01 mK of T C) corresponding to t = ±3 x 10 -7 . 
(b) A maximum of the phase located within 20 /iK of T c . The inset shows that the 
maximum was superposed on a broader minimum. 

17. Model cross section of the CVX electrodes used to calculate electric fields. The field 
near the screen depended on the potentials Vj and Vj, the distances Xj, x%, and a, and 
the wire radius R. 
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